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SUMMARY 


The results o£ an. investigation of the effects of grazing 
flow on the steady-state flow resistance and acoustic ■tm-pedanoe 
of seven Feltmetal and three Rigimesh thin porous -faced liners 
are presented. A state-of-the-art review of previous nongrazing 
flow studies is also presented. The steady-state flow resistance 
of the ten specimens were measured using standard fluid mechanical 
experimental techniques. The acoustic impedance was measured us- 
ing the two-microphone method. The principal findings of the 
study are (1) the effects of grazing flow were measured and found 
to be small, (2) small differences were measured between steady- 
state and acoustic resistance, (3) a semi-empirical model was 
derived that correlated the steady-state resistance data of the 
seven Feltmetal liners and the face-sheet reactance of both the 
Feltmetal and Rigimesh liners. 

These findings suggest that nongrazing flow steady-state 
flow resistance tests could provide considerable insight into the 
acoustic behavior of porous liners in a grazing flow environment. 

1. INTRODUCTION 

Sound absorbent treatments consisting of cavity-backed por- 
ous faced liners are currently being used to control internally 
generated machinery noise. For application to the inlet and ex- 
haust ducting of jet engines, the efficient design of the liner 
should incorporate the effects of both high grazing flow veloci- 
ties and intense sound pressure levels. Unfortunately, there 
have been few investigations of the effects of grazing flow and 
sound pressure on the relationship between the liner impedance 
and its construction (i.e., sheet thickness, porosity, fiber 
diameter, etc) . The purpose of this report is to describe the 
results of a fundamental investigation of the effects of grazing 
flow and sound pressure on the steady-state flow resistance (here- 
in referred to for brevity as d.c. resistance) and the specific 
acoustic impedance (i.e., resistance and reactance) of ten cavity- 
backed porous -faced specimens. The acoustic measurements were 
made using the two-microphone method. This method was used by 
Hersh and Walker^ to study the effects of grazing flow and sound 
pressure level on the impedance of isolated orifices. 

Feder and Dean^ studied experimentally the effects of graz- 
ing flow and sound pressure on the d.c. resistance and acoustic 
impedance of a variety of acoustic liner configurations including 
porous-faced liners. The acoustic measurements were conducted 
using an impedance tube. They concluded that the d.c. and a.c. 
resistance increased and the reactance decreased with increasing 
grazing flow velocity and sound pressure level. The correlation 
between the d.c. and a.c. resistance for three porous-faced mat- 
erials at a fixed grazing flow speed and sound pressure level was 
found to be generally poor; the a.c. resistance values were con- 
sistantly higher (e.g., an a.c. value of 600 ray Is vs. a d.c. value 



of about 370Rayls for a Feltmetal specimen) than the correspond- 
ing d.c. values. Feder and Dean did not interpret their results 
physically. 

P. Dean^ used the two-microphone method to measure the 
effects of grazing flow on the impedance of cavity-backed per- 
forates and porous faced liners. His results were in contrast 
with Feder and Dean's results in that the acoustic resistance of 
some materials did not change while other decreased with increas- 
ing grazing flow; the reactance increased with increasing grazing 
flow. Dean offered no explanation. Plumblee, et. al'' also in- 
vestigated experimentally the effects of grazing flow on porous- 
face materials using the two-microphone method. Since the graz- 
ing flow speeds never exceeded 15 m/sec., they found no signifi- 
cant change in impedance. None of these studies investigated the 
nonlinear behavior of porous-faced materials. 

The present understanding of the connection between porous 
material properties and a.c. and d.c. impedance for nongrazing 
flow is reviewed in Section 2. The results of an experimental 
and analytical investigation of the d.c. resistance and the a.c. 
impedance of ten porous-faced liners are described in Section 3 
and 4 respectively. The major findings of the report are summar- 
ized in Section 5. 
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DEFINITION OF SYMBOLS 
Definition 

cross-section area, (meter^) 

average spacing between fissures (meters) ; 
also constant 

average slit height of fissure (meters) ; 
also constant 

speed of sound (meters/second) 
skin friction coefficient 
numerical coefficient (see Eq . 9) 
drag coefficient 

diameter of porous specimen (meters) ; also 
drag (Newtons) 

frequency (Hz) 

gravitational constant (meters/second’s^) 

height of duct containing specimen (meters) 

complex wave number within porous material 
(meters"^ ) 
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Symbol 

k 

K 

L 

^ch 

m 


M 

N 

Ps 


P'j’ 

P 

P 

P’ 

Ap 

Q 

R 

P-ac 

P-dc 

Re 

r 

S 

t 

T 

u 

V. 

V 


Definition 

material permeability constant (meters^) 

material hydraulic conductivity constant 
(meter/second) ; also compressibility mod- 
ulus (= 3p'/3p'-l/Newtons/meters^) 

'o 

thickness of porous specimen (meters) ; 
also length scale characterizing interac- 
tion between cylinders (meters) ; also 
depth of cavity (meters) 

length scale characterizing thickness of 
unsteady viscous layer (meters) 

material structure factor (ratio of effec- 
tive density of the air in the pore space 
to the actual density of the air) 

material specific internal area (meters"^) 

number of fibers or capillaries 

static pressure upstream of specimen 
(Newtons /meter ^) 

static pressure in test section (Newtons/ 
meter^) 

total pressure in test section (Newtons/ 
meter ^) 

material porosity 

steady-state pressure (Newtons/meters^) 
acoustic pressure (Newtons/meters^) 
pressure drop across specimen (N/m^) 
volume flow, (meters Vsecond) 
average hydraulic radius (meters) 
acoustic resistance (Rayls-mks) 
steady-state flow resistance (Rayls-mks) 
Reynolds number 
radial distance (meters) 
spacing between cylinders (meters) 
time (seconds) 
material tortuosity 

longitudinal component of velocity (meters/ 
second) 

test section speed (meters/second) 
volume of cavity (meters^) 
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Symbol 

Vd 

Vx’ 

X 

X 

y 

w 


e 

a 

a(6) 

3 

6 

r 

p, p’ 

¥ 

Y 

T 

X 

4> 

y 

V 
(jO 
K 

Subscripts 

C )o 

( )c 
inc 

cav 

s 

t 


Definition 

average flow speed incident to porous mat- 
erial (meters/second) 

acoustic velocity (meters/sec) 

reactance (Rayls-mks) 

axial and/or longitudinal coordinate (meters) 

normal and/or lateral coordinate (meters) 

width of test section containing specimen 
(meters) ; characteristic (complex) impedance 
of porous material (Rayls-mks) 

constant defined by Eq, (46) 

constant; also attenuation constant (meters^) 

probability distribution function 

constant; also phase constant (degrees/meters ) 

average cylinder diameter (meters) ; also 
average pore diameter (meters) 

parameter defined by Eq. (26) 

steady-state and acoustic density (kilograms/ 
meters ^) 

velocity potential (meters^/second) 

ratio of specific heat of air 

material thickness (meter) 

wall shear stress (Newtons/meter^) 

sound wavelength (meters) 

phase change across specimen (degrees) 

fluid coefficient of viscosity (kilograms/ 
second/meter) 

fluid kinematic viscosity (meters^/second) 
radian frequency (Hz) 

effective compressibility of the air within 
the porous material (Newtons/meters^) 


denotes open quantity 
denotes closed quantity 
incident 
cavity 

along stream tube direction 
face- sheet 
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Superscripts 

C )* 


Definition 


denotes acoustic quantity 
2. REVIEW OF D.C. AND A.C. IMPEDANCE PREDICTION MODELS 

In contrast to the sparse number of investigations of the 
effects of grazing flow on the d.c. and a.c. impedance of porous- 
faced materials, there have been a large number of investigations 
of their nongrazing flow behavior. A review of the more impor- 
tant papers is given below to provide background information 
necessary to place in perspective the results of the present in- 
vestigation. The d.c. models are described in Section 2.1 
followed by the a.c. models in Section 2.2. 

2.1 D.C. Models 


The review of the d.c. resistance models is divided into 
linear and nonlinear parts. The linear models represent varia- 
tions and/or refinements of Darcy's Law®. Darcy was the first 
to establish a linear connection between the deriving pressure 
gradient Ap across a porous specimen and the average flow speed 
Vp approaching the specimen. For sufficiently high average speeds, 
the connection between Ap and Vp is nonlinear due to a quadratic 
relation between Ap and Vp. Most of the details of the various 
models described below are based on the excellent review by Bear® . 

2.1.1 Linear Models 


In 1856 Darcy investigated the flow of water in vertical 
homogeneous sand filters in connection with the fountains of the 
city of Dijon, France. From his experiments, he concluded that 
the volume rate of flow Q through the filter was proportional to 
the product of the cross-sectional area of the filter A and the 
pressure difference across it (Ap) and inversely proportional to 
the filter length (L) leading to the relationship 



where the constant K is called the material hydraulic conductiv- 
ity, p is the density of the water and g is the gravitational 
constant. Later researchers separated the influence of the ma- 
terial from that of the fluid by defining K = pgk/y where y is 
the viscosity of the fluid and k is the specific permeability of 
the porous material. Substituting this expression for K into 
Eq. (1) yields the modern form of Darcy's Law, 



where Vj) is the average flow speed approaching the porous mater- 
ial . 


A particularly attractive and straight-forward application 
of Darcy's Law is to model a porous material as a collection of 
straight capillary tubes. Applying Hagen-Poisseuille ' s law to 
the flow through a capillary tube of radius 6, Scheidegger^ de- 
rived the following connection between Vp and Ap, 



(3) 


Comparing Eqs. (2) and (3), the permeability constant k = 6^/32 
for a capillary tube. If there are N tubes per unit cross- 
section of the model (normal to the direction of flow) then the 
porosity P = N(7 t 6^/4) and the filter speed through the porous 
material is 



What is important in Eq. (4) is that the filter speed is propor- 
tional to the product of the square of the average tube diameter 
6^ and the material porosity P; the coefficient 1/32 is obviously 
unimportant . 

Scheidegger generalized the capillary model described above 
by replacing the porous material with a distribution of capillar- 
ies of different diameters. He introduced the quantity a(6)d6 to 
represent the fraction of capillaries having radii between 6 and 
6+d6. Assuming that the volume occupied by these capillaries 
parallel, say, to the xi direction, is l/3PAxia (6) dS per unit 
cross-section, the average filter speed Vd . is 



where Vp^. (6) is the average speed in the capillaries of diameter 
6 given by Eq. (3). Comparing Eqs. (5) and (2) yields for the 
permeability constant 


k 


P 

96 



a(<5)d<5 


o 


( 6 ) 


6 



Scheidegger also generalized the above model by permitting the 
capillary tubes to be tortuous (i.e., non-straight) and of var- 
iable diameter. He proposed the following expression for k. 


k=3 




C7) 


Here T is the "tortuosity" factor defined as the ratio of the 
length of the actual flow to the thickness of the porous material. 


Irmay® applied the capillary model to predict 
of narrow fissures (i.e., slits of constant width), 
rock being the closest application of his approach, 
the following expression for the permeability. 


the resistance 
a fractured 
Irmay derives 


k = P®/ 12 (l -P)'' = Pbyi2 


( 8 ) 


where b is the average slit height of the fissure and a is the 
average spacing between fissures. Here the material porosity 
P = b/(a+b) . 

In 1927, Kozeny® applied the concept of hydraulic radius R 
to connect the permeability constant k to the material properties. 
Defining R as the ratio of volume of a conduit to its wetted sur- 
face area, he derived R = P/M, where M is the specific surface 
area of the porous material. Using this definition, R may be 
interpreted as the average hydraulic radius of quite complicated 
flow channels. By solving the Navier-Stokes equations for all 
channels passing through a cross-section normal to the flow, 

Kozeny derived the following expression for the material permea- 
bility constant k. 


k=CoP"My^(l-p)* 


(9) 


where Cq is a numerical coefficient called Kozeny’ s constant 
(Cq = 0.5, 0.562, 0.597, 0.667 for a circle, square, equilateral 
triangle, and rectangle respectively) . 

Recently theoretical estimates of the permeability constant 
have been derived using a drag model. The idea here is that the 
pressure drop through a material is related to internal shear 
(tangential) and normal (pressure) stresses by assuming that a 
porous material can be modelled as a collection of closely packed 
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spheres. For example, Rumer and Drinker^® and later Rumer“ 
derived the following expression for the material permeability 
constant , 


k = 


X(1-P) V L / 


( 10 ) 


where d is the porous material average diameter, X is a coeffi- 
cient that takes into account the effects of neighboring spheres 
(for a single sphere in an infinite fluid X = 3 tt) and 3 is a shape 
factor (3 = tt/ 6 for a sphere). 

Iberall^^ assumed that highly porous materials could be 
modelled as a collection of randomly orientated circular cylindri- 
cal fibers of the same diameter. He then derived an expression 
for the material permeability by assuming that the pressure drop 
necessary to overcome the viscous drag is linearly additive for 
all of the fibers. The idea here is that the average speed is 
sufficiently low, the pressure drop across the material is equal 
to the viscous drag of the fibers. Since Iberall ignored inter- 
actions between fibers, his analysis is restricted to large se- 
paration distances between fibers and hence to high porosity 
materials . 

Iberall assumed that given n fibers of unit length per 
unit volume, n/3 of them will be orientated perpendicular to the 
local upstream velocity Vp. Using Lamb's^^ calculation of the 
drag per unit length for a cylinder perpendicular to a stream, 
Iberall derived the following expression for the pressure drop 
across a specimen 

dp _ (1-P) . 4 - In (Re) . pe = Vn**/ Pv ri n 

“T 3pF 2 -in (Re) ’ °/ 


where 6 is the average radius of the fibers. Comparing Eqs. (11) 
and (21) the permeability k is 


3P5^ 2- In (Re) 

16(1-P) 4 -In (Re) 


( 12 ) 


Equation (12) is unusual because it shows that the material per- 
meability is not only a function of the material properties but 
is also a function of the fluid inertia through the Reynolds 
number. Scheidegger^ suggests that the drag theory may represent 
a generalization of Darcy's Law. 
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2.1.2 Nonlinear Models 


For the linear relationship between Ap and Vp to be valid, 
the Reynolds number characterizing the flow through the porous 
material must be extremely small. As used here, the Reynolds 
number is defined as Re = Vj)6/Pv where Vp/P is the average velocity 
in the material and 6 the average pore diameter of the material. 

The Reynolds number can be interpreted as a measure of the rela- 
tive importance between inertia and viscous forces. To show 
this, consider flow over an isolated cylinder of radius 6. By 
nondimensionalizing all length changes by the cylinder diameter 
6 and all velocity changes by the toeai upstream velocity Vp/P, 
it follows that the order-of -magnitude of a typical nonlinear 
term pu8u/3x is Cp/5)(Vp/P)^ and that the order-of-magnitude 
of a typical shear stress term is /j3^u/8y^ ~ pVp^/vP^. The ratio 
of the nonlinear inertial to viscous terms is Vp6/Pv=Re. 

It is clear that the flow within a porous material cannot 
be accurately modelled by a collection of straight, constant area, 
capillary tubes. A more realistic model would have to incorpor- 
ate both twisting fully three-dimensional motion with the fluid 
accelerating and decelerating within the material. Recognizing 
this Scheidegger believes that the onset of nonlinearity first 
occurs due to laminar flow turning and fluid accelerations and 
decelerations. At higher speeds, nonlinear losses arise due to 
turhulenoe generated within the porous material. According to 
the review by Bear, initial deviations from Darcy's linear law 
occur for Reynolds number as low as Re = 0.1. The deviations 
increase for increasing Reynolds number. These deviations are 
believed to arise from the onset of laminar nonlinearity as dis- 
cussed above. The onset of turbulence generated nonlinearity 
does not occur until Reynolds numbers between 60 to 150. In 
general, the following relationship between Ap and Vp character- 
ize the d.c. resistance. 


Ap = a Vp + b Vp^ ; a, b constants 


(13) 


Here a is proportional to y and independent of the fluid density 
p and b is proportional to p and independent of y. 

Blick^'' analyzed the nonlinear behavior of porous materials 
in a rather interesting way. His model consisted of a bundle of 
capillary tubes with orifice plates spaced along each tube, sep- 
arated by a distance equal to the mean pore diameter. Assuming 
the fluid is homogeneous, Newtonian and one-dimensional, Blick 
derives the following expression relating p and Vp. 
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2 


(14) 


liE. = 

2cf Re|/ 

Vn 4- 


ax 

P(5" 

D 4- 

2bP^ 


where C£ is the skin friction coefficient of the flow along the 
tube walls of constant diameter 6, Re=Vj)6/vP is the diameter 
based Reynolds number and is the drag coefficient of the ori- 
fice plate. Blick's model is important because it shows quali- 
tatively the nature and form of nonlinear pressure losses. For 
a detailed discussion of other models, the reader is referred 
to the review by Bear. 

This review of flow nonlinearity concludes by quoting from 
Bear's review. "Cheaveteau and Thirri-ot^^ perform experiments in 
several two-dimensional models that permit visual ovservations of 
streamlines and streaklines . In a typical models which has the 
shape of a conduit that diverges and converges alternately ^ they 
observe streaklines. For Re<2, the flow obeys Darcy's law and 
the streamlines remain fixed. As Re increases ^ streamlines start 
to shift and fixed eddies begin to appear in the diverging areas 
of the model. They become larger as Re increases . At Re=75 
turbulence appears and starts to spread out as Re increases . 
Turbulence covers some 50% of the flow domain at Re=115 and 100% 
of it at Re=180 . The deviation from Darcy's law is observed at 
Re=2-3. Thus the deviation from Darcy's law as the velocity in- 
creases is associated with gradual shifting of streamlines due 
to the curvature of the microscopic solid walls of the pore 
space . " 

2.2 Acoustic Models 


Despite widespread application of porous materials as sound 
absorbers, surprisingly little fundamental research has been con- 
ducted in the past twenty years to understand their acoustic be- 
havior. Most of the recent research has been applications orien- 
tated in the sense that existing parameters used to define the 
important physical properties, namely structure factor, bulk mod- 
ulus of compressibility and dynamic flow resistance are measured 
or assumed known. The results of three investigations are re- 
viewed below representing the pioneering work of Scott Zwikker 
and Kosten^^ and Morse^®. These three are believed to adequately 
summarize the state-of-the-art. None of the investigators consid- 
ered the acoustic nonlinear behavior of porous materials. 

Scott was one of the first to be concerned with the propa- 
gation of acoustic disturbances in porous materials. He assumed 
the porous material to be homogeneous, isotropic and sufficiently 
large in extent so that reflected waves were negligible. He 
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further assumed that the porous material consisted of very many 
small fibers or particles packed together so as to leave between 
them interconnected air-spaces or pores of irregular shape. De- 
fined in this way, the linearized equations of motion for sound 
waves propagating through the pore space are obtained. The lin- 
earized conservation of momentum equation becomes 

-V7p' =mp„ + Ra.c. 


where m is the ratio of the effective density of the air in the 
pore space to the actual density Po of the air (called the struc- 
ture factor by other investigators), p' is the acoustic pressure, 
is a dynamic resistance constant analogous to the d.c. 
resistance used in Darcy's law, P is the material porosity and 
Vp is the vector velocity approaching the porous material. The 
corresponding linearized mass conservation equation is given by 

''•( Yp )=0 (161 


By introducing the acoustic velocity potential Vp=grad ijj , Scott 
derives the following linearized wave equation 


+ h"'P =0 


(17) 


with 


iRcic-Pco 

K K 


(18) 


The quantity h is the complex propagation constant within the mat- 
erial and K is the effective compressibility of the air in the 
pore spaces. By introducing the complex density p' defined as 


P'= 

P(j^ 


(19) 
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Scott derived the wave equation in a form similar to the wave 
equation in free space. The remainder of his paper describes 
the use of a probe tube to measure the real and imaginary parts 
of the propagation constant h. Writing h=3-ia, the decay of 
intensity with distance provides a measurement of the attenua- 
tion constant a. The corresponding change in phase of the 
acoustic pressure with distance provides a measuremnt of the 
wavelength, and therefore of 3. With h specified, the complex 
density p' follows immediately from Eq. (19)* 

Scott compares his theory with measurements of a variety 
of rock-wool porous samples of different thicknesses. At 
f=200 Hz, he concludes that for a rock wool sample, the effec- 
tive density of the air trapped within the pores was roughly 
twice its free air value; it decreased to about 1.32 at 
f=4000 Hz. His measurements also showed that at f=200 Hz, the 
dynamic resistance was equal to its d.c. value. The speed of 
sound varied from Cq/y at low frequencies to Cq at the higher 
frequencies, thus indicating that low frequency sound propagat- 
es isothermally and high frequency adiabat ically . Scott's paper 
is obviously important. His desire to recast the sound propaga- 
tion equations in a form equivalent to that in free space is 
understandable. Unfortunately, in the process of forcing the 
analogy, he as well as Zwikker and Kosten may have induced work- 
ers in the field into believing that the effective mean density 
of sound waves propagating in porous materials is a function of 
the structure factor m, porosity P, flow resistance and 

frequency to. In reality the mean density of the air is constant 
since the air can be considered incompressible. While m, P, r 
and to affect the sound pressure, density and velocity, they do 
not affect mean density. This is most evident at low frequencies 
where the a.c. resistance is equal to the d.c, flow resistance. 

At the low particle speeds, generated by the sound waves, the 
incident sound driving pressure cannot increase the air-density - 
clearly the process is incompressible. Thus interpretations 
other than increases in the fluid mean density are required. 

In 1949, Zwikker and Kosten published their famous book 
entitled, Sound Absorbing Materials, which dealt in part with the 
connection between the quantities h and p' and the material struc- 
tural properties. Zwikker and Kosten (herein referred to Z § K) 
introduced the following assumed solution describing the propaga- 
tion of a plane sound wave in a homogeneous isotropic medium ex- 
tending to infinity. 


P 



-ax + i 
Ae 


(wt-^x) 


( 20 ) 
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where a is the attenuation constant o£ the sound wave, 3=to/c is 
the corresponding phase constant (c is the phase velocity of the 
sound wave in the material] and A is an arbitrary amplitude. A 
second constant is defined, called W, the characteristic (complex) 
impedance of the infinite extent material, W=p’ (x, t) /v^’ (x, t) . 

Here Vx*Cx,t) is the sound particle spee d in the x - direction 
within the material. W is defined as /1/pJC \thexe K= (3p * / 3p* ) /Po is 
related to the complex speed of soimd in the material. W is com- 
plex due to internal material damping of the propagating sound 
wave; it denotes that acoustic density variations within the mat- 
erial are not in phase with the driving pressure variation. 

There are two limiting cases in which p’ and p* vibrate in phase 
with other. In the first limit, the sound frequency is suffi- 
ciently low that the transmission of heat to and from the solid 
frame is rapid enough so that the enclosed air is kept at constant 
temperature (isothermal case). At high frequencies, the trans- 
mission of heat is so slow that the air vibrates adiabatically . 
Between these limits hysteresis, occurring between the instantan- 
eous heating and cooling of the air by the sound field and the 
subsequent transfer of heat to and from the material, results in 
a frequency dependent complex propagation constant. 

The continuity and momentum conservation equations govern- 
ing the motion of the sound waves are, respectively. 


3x 


- ^ Z' \ ^ and liE 
p^VdpV at ax 



at 


R Vy' 
’ a.c. X 


( 21 ) 


where P is the material porosity, Rac the material a.c. resistance 
an d m th e structure constant. For low frequency oscillations 
(/w6^/v<<l) Ra.c.=^^d.c. quantity Rd.c. follows from Poise- 

uille's law governing the behavior of steady- state flow in capil- 
lary tubes, namely 


Rdc = p5' 


( 22 ) 


where 6 is the avera ge material pore radius. For high frequency 
oscillations (/o)6^/v>>l) , 


R,. = 


m 

P 




(23) 
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Z § K offer the following interpretation of the term m/P 
that multiplies the density Pq of the free air in the momentum 
equation. They split the factor into two parts, m and 1/P. 

The quantity m is thought to consist of two parts, one of which 
is related to the internal friction of the air, the other to the 
orientation of the material pores. A third part also exists, 
the inertia of the structural material, but Z § K feel that in 
most cases this effect is small. 

Of the various contributors to m the pore orientation is 
the most important. Z § K show that it contains the factor 
(l/cos^6). For randomly orientated pores m^3 . According to ex- 
periments, m has values lying between 3 and 7. By introducing 
the apparent (complex) density defined as 


f 

P = 


O *“ ^ ^ d.c. 

P 0> 


(24) 


the equations of motion for a traveling wave is brought into the 
universal form. 


■T p = icjp'Vx 


(25) 


where 


P = \u^py K , W = ^Kp' 


(26) 


Zwikker and Kosten's book is an accurate state-of-the-art 
summary of the acoustical understanding of porous materials 
through 1949. Following Scott, the authors chose to cast the form 
of the equations of motion into a universal form - i.e., similar 
in form to the equations of motion describing the propagation of 
sound in air. Thus their interpretations suffer the same limita- 
tions as those of Scott. 

Morse extends the approach adopted by Zwikker § Kosten 
(first presented by Rayleigh) and by Scott, Morse and Bolt^“ , 
Beranek^^ § others to model the propagation of sound in granular 
media. He assumes that all of the grains have roughly the same 
diameter. He derives the following equations governing the sound 
motion in the granular material. 
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V- u'= -p IE. 

at 


( 27 ) 


P’ 


^ 2 _ I 

Co P 




at 


_ R 


V 


Here v' is volume o£ fluid per unit area in direction pressure 
gradient 

p ' , Po are the acoustic and equilibrium fluid density respective- 
ly 

p' is the excess instantaneous pressure 
m is the structure factor (explained below) 

P is the material porosity 
Rdc is the dynamic flow resistance 


Morse defines the structure factor m as arising because of 
the way in which the solid material constrains the motion of the 
fluid. Not all the fluid is 'in the direction of the pressure 
gradient VP'. To explain this, consider small tubes inclined at 
an angle 0 to the pressure gradient. 



is 


In the S-direction (along the tube) , the momentum equation 


__3p' 
3t as 


(28) 


Now Us'cos0=ux' 
yields 


and scos0=x^3/3s=cos03/3x. Thus substitution 
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3vx _ - 3p' 
9t 3x 


Averaging the tubes over all possible angles assuming that all 
the tubes are orientated uniformly in direction, the probability 
that a given tube lies between the spherical angles 9, Q+d0 and 
(p, <p+d(p is unity, then 


<cos^6> 


1 

27T 



1 • cos* 9 • sin 6 d d 



o o 


Thus, m=3. Morse does not allow for inertial coupling between 
the motion of the porous material and its effect on the air. 

Assuming a plane wa ve solution of the form p ' (x ,t) 

Morse shows that k= (w/cq) ’^m- iRdcP/Po^^* Assuming further that 
RdcP/Po"^^l> "the propagation constant k may be written 


k 




2Po Cc 


(30) 


Thus the absorption (a) and sound propagation speed (c) 
within the material follows immediately to be 




« - Rd.c.P/ 2p„c/ vm 

c = ! 


(31) 


Clearly, the structure factor m(m>l) as defined above results in 
a reduced sound propagation speed. Obviously m includes heat 
transfer effects (isothermal, adiabatic) not explicitly discussed 
by Morse. 


Morse proceeds to connect the dynamic resistance R<jc to the 
known acoustic behavior in small tubes . At low frequencies where 
the acoustic Reynolds number Ra = /u)6^ /v<<l , Hj^,-Rg^^=8vm/P6 ^ . In- 
serting R(Jq into E g. (31) shows that a~v^ At high frequencies 
where Ra>>l, Rac^ni/ZytopQ/PS . The absorption a is also proportion- 
al to Morse concentrates only on the high f requency case. 

He analyses available data to show that a~/vR^. 
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Morse's paper is attractive because he does not identify 
the quantity m/P as increasing the effective density of the air 
trapped within the porous material as do Scott and Z § K. In- 
stead, he interprets m to affect the magnitude of the pressure 
gradient because of lateral fluid motion induced by the structure 
not in the direction of the pressure gradient. The porosity P 
is obvious and requires no explanation. Morse interprets the 
effective speed of sound c in the material to be c=Co/»'^where 
Co is its adiabatic value in free air (recall that Morse consi- 
ders only the high frequency solution) . While this interpreta- 
tion is valid at high frequencies, it is incorrect at low fre- 
quencies wherein Cq should be redefined as co/y (the isothermal 
sound speed) . 

3. STEADY- STATE FLOW RESISTANCE TESTS 

Many investigators have shovm that in the absence of graz- 
ing flow, the d.c. and a.c. resistances of porous materials are 
equal (e.g,, see the study by Zorumski and Parrott^^). For 
application to the aircraft industry, it is of interest to deter- 
mine whether this equivalency is still valid in the presence of 
a grazing flow. Accordingly, d.c. flow resistance tests were 
conducted with the porous-faced liner specimens. 

Table I summarizes the manufacturers* predicted acoustic 
properties. Specimens #1-7 were made from Feltmetal fiber metals 
and specimens #8-10 from Rigimesh stainless steel wire cloth. 
Photographs enlarged to illustrate the detailed structure of 
typical Feltmetal and Rigimesh specimens are shown in Fig. 1. 

3 . 1 Experimental Apparatus 

The effects of grazing flow were studied experimentally in 
the HAE wind tunnel shown in Figure 2. The test section, made 
from .0127 meters thick transparent acrylic sheet, has a remov- 
able side-wall containing the test specimen. The porous test 
material was backed by a cylindrical tube of diameter .0508 
meters; it was secured and flush mounted to the test section side 
wall by means of three set screws. Flow through the side branch 
was controlled by means of blowers used in the suction mode 
(i.e., from the test section into the side branch) or blowing 
mode (i.e. , from the side branch into the test section). The 
side branch volume flow rate was monitored by means of sharp- 
edge orifices. By accurately measuring the pressure drop across 
the orifices with Charles Merriam Company calibrated slant mano- 
meters, volume flow rates as low as 8.5 x 10’® cubic meters per 
second were recorded. 

Test section speeds up to 91.44 meters/second were measured. 
The velocity profile for V^=85.6 meters/second, is shown in Figure 
3. As shown, the boundary layer is turbulent - measurements 
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showed it to be turbulent for all test speeds considered. For 
distances less than .00051 meters from the tunnel wall, the data 
indicates that the flow is separated due to interaction between 
the boundary layer probe and wall. The data matches closely the 
classical l/7th velocity law profile for .00051 meters<y<. 00762 
meters . 

The d.c. resistance R of a typical specimen mounted as 
shown in the lower part of Figure 2 is defined as the ratio of 
the pressure drop Ap across the material to the average (within 
the .0508 meter diameter tube) velocity Vp, 


R 



The pressure drop was measured directly from a calibrated Charles 
Merriam slant manometer. The average speed was calculated using 
standard ASME techniques'^ from the pressure drop measurements 
across the orifices (orifice diameters of .00635 meters, 0127 
meters and .01905 meters were used to extend the useful range of Vp 
from a minimum of .04 to a maximum of 10 meters/secord) . 

3 . 2 Test Results 

Figures 4-9 summarize the effects of grazing flow on the 
d.c. flow resistance of the ten specimens operating in both the 
suction and blowing mode. In agreement with the discussion pro- 
ceeding Eq. (13) , the data can be accurately presented by the 
empirical curve fit 


R 


d.G. 



a + b Vp 


(32) 


The various measured values of a and b are tabulated in Table II. 
Inspection of Table II shows that, in general, there are only 
small differences between the blowing and suction data. Further, 
the effects of grazing flow are also shown to be small. Most of 
the differences are of the order of 201 or less. These differ- 
ences are considered small relative to the effects of grazing flow 
on the d.c. resistances of orifices as measured by Feder and Dean. 

The measured negligible differences in flow resistance for 
the specimens operating in the blowing and suction modes indicate 
that the material permeability is isotropic in the axial direction. 
Flow nonlinearity is shown to become important for Vp>2 meters/ 
second for specimens 2, 3, 4, 8, 9 and 10 for Vp>.30 meters/second 
for specimens 1, 5, 6 and 7 (recall from Table I that these 
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specimens should be equivalent) . This is reflected in the 
various values of the coefficient b summarized in Table ‘II. 

Since b is a measure of the importance of flow nonlinearity, 
following the arguments of Scheidigger, it is a function of the 
area changes associated with the various pore diameters. The 
idea here is that relatively large pressure gradients are re- 
quired to balance the relatively large convective velocity 
gradients which arise from abrupt changes in internal pore 
diameter cross-sectional areas. 

Specimen 1, 5, 6 and 7 should, in principle, exhibit 
roughly the same steady-state flow resistance behavior. Clear- 
ly , the data characterizing specimens 5 and 6 are reasonably 
close - the coefficients a and b are within 20%. The corres- 
ponding coefficients characterizing specimens 1 and 7, however, 
are quite different. These differences are probably associated 
with the variability inherent with their manufacture. The de- 
tails of this variability are not well understood nor documented. 

The negligible changes in flow resistance with grazing 
flow require special comment. A small positive pressure drop 
APs=Ps~P^ induced across the specimen by the grazing flow was 
observed. The flow control valves across the side branch were 
closed so that there was zero mean flow across the specimen. 
According to the investigations by first Beavers and Joseph^"* 
and later Taylor^®, the observed pressure drop arises because of 
a local recirculating (vortical) flow induced across the speci- 
men. A sketch showing the vortical flow pattern is shown in 
Fig. 10. The recirculating flow is driven and maintained by 
steady-state shear stresses transferred across the porous sur- 
face. Figure 11 shows the effect of correcting the data by 
subtracting this pressure drop from the measured blowing data 
(adding for the suction data) for the data of Specimen #4. With- 
out this correction, the flow resistance would increase without 
limit as Vj)->0 as indicated in Fig. 11 (in the suction mode, it 
would decrease without limit). Table III summarizes the measured 
pressure drops APs in meters of water across the ten specimens 
at a grazing flow speed of V^=79.6 meters/ second . The measure- 
ments showed that (APs) is linearly proportional to the test sec- 
tion head (Py-P^) where Py is the test section total pressure and 
P^ is the test section static pressure. 

3 . 3 Derivation of Semi-Empirical Prediction Model 


A simplified prediction model of the nongrazing linear d.c. 
resistance of thin porous-faced liners is derived. The intent 
here is to correlate the linear d.c. resistance data. It is under 
stood because of the variations inherent in the manufacturing pro 
cess that only an approximate correlation is possible. The photo 
graphs of the Feltmetal specimens, shown in Figure 1, suggest that 
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they can be modelled as a collection o£ small diameter cylinders. 
A model o£ the Rigimesh specimens will not be derived because 
o£ the very limited number o£ specimens supplied by the manu£ac- 
turer. To simpli£y the model, the £ollowing assumptions are 
made. The cylinders have mean diameter 6. They are aligned 
perpendicular to the incident £low Vp as shown in Fig. 1. The 
specimen has zero thickness, that is, all the cylinders are as- 
sumed to lie on a plane surface - this treats the specimen sur- 
face as a momentum sink in a manner similar to that proposed by 
Zorumski and Parrott. The cylinders are separated by an average 
center- to -center spacing S as shown schematically in Figure 1. 
Here the specimen is placed in a duct of width 
with incident flow Vp as shown. The parameters 
related as follows 


H NS for N »1 


Consistent with the liner zero thickness 
porosity P is defined as the ratio of the open 
total area Ay, 

p= Aq - 1 _ NdWy^HW 
At 


W and height H 
H, N and 6 are 


(33) 


assumption, the 
area Aq to the 


(34) 


Substituting Eq. (33) into Eq. (34) and assuming N large yields 
the following approximate estimate of the porosity. 


P 



(35) 


The d.c. resistance is calculated by assuming that the cy- 
linders act as a momentum sink. From momentum balance considera- 
tions, the following relationship exists between the pressure 
drop Ap across the specimen and the drag force due to the cylin- 
ders , 


dp • H W = ND 


(36) 


where N is the total number of cylinders (of width W) and D is the 
drag force on each cylinder. Now assume for very slow upstream 
average speeds Vp, that the cylinder drag is proportional to the 
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skin friction shear stresses where 


D = 


TTdW- 


(^p/ p) 


"ch 


TT^W 


C37) 


Here C^/P) is the local pore velocity incident to the cylinder and 
Lch is a measure of the thickness of the viscous shear layer 
normal to the cylinder surface. For very low Reynolds number 
flows, the characteristic length Lch for an isolated cylinder is 
the diameter 6. For collections of closely spaced cylinders, how- 
ever, L^h niust be a function of the separation parameter S. Thus 
for S/6>>1, Lch=<5 and for S/6<1, Lch“f (S/<S) • Here S/6 may be 
thought of as an interaction parameter. The following simplified 
empirical relationship is suggested. 


as 


"ch- 


1 + 


a (s/^ 


C38:) 


For S/6>>1, Lch~^- P°r S/6<<1, L^h-aS where a is a constant that 
must be determined from experimental data. Hopefully, a will be 
independent of specimen geometry. Combining Eq. (37) and Eq. (38) 
and substituting into Eq. (35) yields 

±12. ~ NTTjt/^pgWP 
HW HWPotS 


1 + 


a (s/ «) 


C39) 


Substituting Eqs. (33) and (34) into Eq. (39) and solving for the 
ratio Ap/Vp yields 

ad \ P / 


The following form of the d.c. flow resistance model is written as 

(7) 


1 - P-i- 


a 


(40) 


where a and 8 are unknown constants. 

The results of applying Eq. (41) to determine the constants 
a and 8 are summarized in the fifth column of Table IV. The best 
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fit to the data occurred with a=0.8; the average value of 
86 4 . 

The model correlates the Feltmetal data quite well over a 
wide range of porosities (0.45 to 0.97), cylinder diameters 
(4x10"® to lO"** meters) and d'.c. resistances (115 to 890 
Rayls-mks) . No particular significance is attached to the values 
of a and 3. The porosity values shown in Table IV were estimated 
by computing the relative densities of the specimens. This in- 
volved subtracting out the weight of the screens and weighing 
the specimens and measuring their volume. The agreement between 
manufacturers estimated porosities and measured porosities are 
quite close. 

No attempt was made to model the nonlinear behavior of the 
specimens. The reason for this is that specimens #2, 3 and 4 
were just barely entering the nonlinear regime at the higher test 
filter speeds. Based on the flow visualization tests performed 
by Chauveteau and Thirriot (see Section 2.1), nonlinear devia- 
tions from Darcy's law were observed at values of Reynolds number 
Re=V]^6/Pv between 2 and 3. Using these values of Re, the last 
column in Table IV shows the values of filter speed Vp correspond- 
ing to Re=2 and 3. Comparing these values with the observed 
values of Vp representing the onset of nonlinearity (the sixth 
column in Table IV) indicate fairly good agreement with Chauve- 
teau and Thirrot's criteria. Accepting Chauveteau and Thirriot 's 
criteria, it is clear that the onset of nonlinearity can be ex- 
tended to very high average speeds Vp by using porous materials 
made of very small fiber diameters. 

4. ACOUSTIC IMPEPANCE TESTS 

In the previous section, the effects of grazing flow on the 
d.c. resistance of ten porous-faced materials were measured and 
found to be negligible. In this Section, the two-microphone 
method is used to measure the effects of grazing flow, sound pres- 
sure level and frequency on the acoustic resistance and reactance 
of the same ten specimens. The test apparatus and the two-micro- 
phone method is described in Section 4.1. The test results are 
summarized in Section 4.2. A non-grazing flow semi-empirical 
impedance prediction is derived in Section 4.3. The model is 
used to correlate the nongrazing test data of Section 4.2. 

4 . 1 Two-Microphone Method and Test Set-up 


The acoustic resistance, total reactance, and face- sheet 
reactance of the test specimens were measured using the two mic- 
rophone method described by Dean®. The method requires the 
specimens to be cavity-backed. Figure 12 shows a schematic of 
the test set-up and the instrumentation used. The acoustic 
resistance and reactance of cavity-backed specimens can be 
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described by the following expressions 
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where SPL(inc) -SPL(cav) represents the sound pressure level dif- 
ference (in dB) between the incident and cavity sound fields and 
(j) represents the corresponding phase difference. The quanties (o, 
L and c represent the sound radian frequency, cavity depth and 
fluid sound speed respectively. The above measurements are ob- 
tained by flush mounting one microphone at the cavity base and 
the other flush with the side walls containing the orifice as 
shown in Figure 12. It is important to locate the incident mic- 
rophone sufficiently far from the specimen to avoid near field 
effects. The microphone should be located sufficiently close, 
however, so that the separation distance is small relative to the 
incident sound wavelength. This is necessary to insure accurate 
measurements of the incident sound wave amplitude and phase. A 
discussion of the errors associated with using the two-microphone 
method is contained in Appendix A. 

A schematic of the instrumentation used to conduct the ex- 
periments is shown in Figure 13. To generate incident sound 
pressure levels up to 150 dB, a JBL type 2480 driver capable of 
producing in excess of 10 watts of relatively "clean” acoustic 
power is used as the sound source. The .0508 meter diameter 
driver throat is coupled to the test section by means of a .0508 
meter to .1016 meters diameter exponential expansion, JBL type 
H-93. Sound pressure levels in excess of 150 dB exceed the input 
capability of the GR-1560-P42 preamp. A 10 dB microphone attenu- 
ator, GR Type 1962-3200 has been added, which extends the 
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measurement range accordingly. 

The signal generated by the Heath lG-18 audio generator is 
amplified by the McIntosh MC2100 100 watt/channel power amplifier 
to power the JBL driver. The audio generator provides a tracking 
signal for the AD-YU Synchronous Filter and phase meter system. 

The 1036 system filters the two microphone input signals to the 
tracking signal frequency with a bandwidth of 5 Hz. The AD-YU 
Type 524A4 Phase Meter reads phase angle between the signals 
independent of signal amplitudes. The phase angle output is 
displayed on the AD-YU Type 2001 digital volt meter. A General 
Radio-1565 1/10 octave filter together with a Heath Type 1M2202 
DA/M is used to record the output signals from each of the two 
microphones. Also the two signals are observed on a TEK 533 
Oscilloscope to visually note approximate phase and distortion 
effects . 

The output of the incident microphone channel of the syn- 
chronous filter is used as a control voltage for an automatic 
level control amplifier. This control amplifier adjusts the drive 
level to the power amplifier in such a way as to keep the incident 
level constant, independent of frequency and amplitude response 
irregularities in the loudspeaker and tunnel. 

As a convenience, a triple ganged 5 dB per step ladder 
attenuator is used to simultaneously increase the power amplifier 
drive level and decrease the synchronous filter input signals so 
that the control loop of the automatic level control amplifier 
always has the same gain. This has the added advantage keeping 
the levels at the AD-YU Filter input constant for all testing 
levels. Since the AD-YU Filter displays a small amplitude -phase 
dependency, this automatic level control improves accuracy as well 
as speed of data acquisition. A phase tracking test of both mi- 
crophones mounted flush in the wind tunnel wall indicated accurate 
phase tracking within +0.2° over a sound pressure level range of 
70-150 dB. 

4.2 Test Results 


Figures (14) -(16) summarize the effects of grazing flow and 
sound pressure level on the acoustic resistance of specimens 
#1-10. Grazing flow is shown to increase only negligibly the zero 
grazing flow resistance data. Only grazing flow data correspond- 
ing to incident sound pressure levels greater than 110-120 dB are 
presented to insure an adequate signal- to-noise ratio. In agree- 
ment with the steady-state flow resistance* data, the zero grazing 
flow data show that nonlinearity effects are important only for 
the Fibermetal FM 134 porous material - specimens #1, 5, 6 and 7. 
Nonlinearity effects are unimportant for specimens #2, 3, 4, 8, 9 
and 10 . 
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The effects of grazing flow and sound pressure level on 
the acoustic reactance of specimens #1-10 are summarized in 
Figures (17) -(25). In all cases, the porous face-sheet reac- 
tance and total system reactance (i.e., face-sheet and cavity 
reactance) are either unaffected or changed by only small (but 
measurable-like 10% differences) amounts from their zero graz- 
ing flow, low sound pressure level values. The frequency at 
which the reactance measurements were conducted requires special 
comment. For each specimen tested, the frequency of the inci- 
dent sound field was selected so that at zero grazing flow and 
at an incident sound pressure level of 90 dB, the phase shift 
across the specimen face-sheet was ninety degrees. Under these 
coditions, the total reactance of each specimen was zero corres- 
ponding to resonance. 

Figures (26) -(31) summarize the acoustic response of speci- 
mens #1, 2 and 7 to variations of frequency. In general, meas- 
urements show that the acoustic resistance is independent of 
frequency as shown, for example, in Figure 26. Specimens #2 and 
7 show, however, an interesting dip in resistance in the neighbor- 
hood of 700-750 Hz. The dip occurs regardless of the amplitude 
of the incident sound or the speed of the grazing flow. Impact 
measurements showed that mechanical resonance corresponding to 
the first resonant clamped mode occurred at 1550 Hz for specimen 
#1, 730 Hz for specimen #2, and 360 Hz for specimen #7. This 
accounts for the observed resistance dips for specimens #2 and 7- 
it corresponds to the excitation of the simple-harmonic fundamen- 
tal mechanical resonance of the porous face-sheet material of 
specimen #2 and possibly to a higher harmonic of specimen #7. 

The dip in resistance physically corresponds to the sound par- 
ticle field converting less mechanical energy into heat within 
the porous material. 

The corresponding variation with frequency of the face- 
sheet and total reactance of specimens #1, 2 and 7 are shown in 
Figures (29) thru (31) . Figure 29 shows that both the face-sheet 
and total reactance of specimen #1 to increase almost linearly 
with frequency. Figures 30 and 31, however, show the face-sheet 
reactance of specimens #2 and 7 to increase abruptly at the fre- 
quency corresponding to mechanical resonance. The total reac- 
tance of these specimens were only slightly perturbed at mechani- 
cal resonance. This is more clearly shown in Figure 32 which 
shows the sound pressure amplitude ratio across the specimen to 
increase significantly for specimen #2 with virtually no increases 
for specimen #1 and 7. Figure 32 also shows that specimen #2 
experiences a large phase change between 700-750 Hz; however, 
specimens #1 and 7 do not. 

To support the nongrazing flow two-microphone test results, 
the impedance of the ten porous samples were measured in a 1.83 
meter long, .0508 meter diameter impedance tube. Since the 
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impedance method is well known it will not be described herein. 
The impedance tube test data are compared with the nongrazing 
flow d.c. resistance and two -microphone impedance test data in 
Table V. Eere the two-miarophone data are eorreoted for aooustio 
radiation from the specimens using the method described in Appen- 
dix A. Comparison of the d.c. and a.c. resistance is generally 
quite good. To first order, d.c. and a.c. resistances are equal. 
Comparison of the reactance data between the two-microphone and 
impedance tube measurement methods is only fair. It is extremely 
difficult to measure accurately using an impedance tube the reac- 
tance of a specimen whose resistance is close to pc. Thus the 
two-microphone method of measuring reactance is believed to be 
superior to the impedance tube method. 

4 . 3 Derivation of Non-Grazing Flow Impedance Prediction Model 

A prediction model is derived of the impedance of thin 
porous-faced cavity-backed liners. Since the liner thickness is 
very small compared to the sound wavelength, the model will be 
derived using the concepts of lumped elements (i.e., the slug- 
mass model). This method was used by Rayleigh^® in deriving the 
impedance of the Helmholtz resonator. Rayleigh’s approach is 
non- fluid mechanical but gives the actual acoustic impedance 
characteristics for low sound pressure levels (i.e., the ’’linear" 
impedance regime) when an empirical end correction is added to 
the slug mass. 

Rayleigh derived the following expression for the impedance 
of a single orifice of diameter d and thickness x backed by a 
cavity of length L, 


Z - Ra,c./P+ i 


p(jJ (t .85 d) - pc cot 



orifice 

resistance 


orifice inertia cavity stiffness 
reactance reactance 


(44) 


where P is defined as the ratio of orifice area to cavity cross- 
sectional area. The impedance defined by Eq. (44) has been re- 
ferenced to the cavity cross-sectional area. The impedance de- 
fined by Eq. (44) can be used to model the impedance of the 
porous materials by replacing the orifice a.c. resistance with 
the d.c. resistance of Eq. (41) with a=0.8 and 3=1865 (the 
data described in Section 4.2 showed this to be a good approxi- 
mation) and by rewriting d=/?D (for the cylindrical cavity shown 
in Figure 12, P=(d/D)^) leads to the following estimate. 
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Here the constant 0.85 in Eq. (44) is replaced with an unknown 
(hopefully) constant e to be determined from experimental data. 
The various values of e computed using Eq. (45) and the two- 
microphone measured values of face-sheet reactance are summar- 
ized in Table VI. With the exception of specimens #5, 6 and 7, 
e is seen to be constant with average value ey^y==0.52 indepen- 
dent of specimen properties for both the Feltmetal and Rigimesh 
specimens. The lack of agreement between specimens #5, 6 and 7 
and the remaining specimens cannot be explained. The excellent 
agreement among the remaining specimens is encouraging and 
suggests that specimens #5, 6 and 7 (which should be equivalent 
to specimen #1) may be improperly fabricated by the manufacturer. 

The final expression for the face- sheet reactance, valid 
for all ten specimens, is 



The substitution for the Feltmetal specimens of the d.c. 
for the a.c. resistance is valid only for low frequencies. A 
physical explanation for the low frequency restriction is des- 
cribed below. 

The steady state model described in Section 3.3 will now 
be extended to unsteady flows. What is required is a reasonable 
estimate of the fluctuating drag of a cylinder immersed in a 
sound pressure field. It is assumed that most of the forces 
exerted by an oscillating sound field on a cylinder occur in the 
form of fluctuating wall shear stresses Thus, the drag force 

per unit cylinder length is balanced by a shear force per unit 
cylinder length given by 



(47) 


Here p is the fluid coefficient of viscosity, 6 is the cylinder 
diameter^ L^h is the extent to which the unsteady shear stresses 
diffuse laterally from the cylinder surface and v^' is the ampli- 
tude of the fluctuating sound particle velocity. 

An order-of-magnitude estimate of the high frequency solu- 
tion to Eq. (47) is described below. The estimate starts with 
the boundary- layer approximations to the momentum equation for 
an unsteady flow past a stationary cylinder. 
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The boundary-layer assumptions require that the lateral region of 
unsteady diffusion of vorticity be small w.r.t the cylinder 
radius. Physically this requires high frequency sound oscilla- 
tions. Linearizing Eq. (48), and solving for the pressure grad- 
ient yields 


I f 2 I 

£l=_p + p ^ 

dx 6t dy^ 



Far from the surface of the cylinder (see sketch) , that is for 
sufficiently large values of y, the particle velocity must satis- 
fy the following boundary conditions. 


Lim 

y — » CO 


dv 

dy 


2L =0 and 


Vy, - ve 


ia>t 


(50) 


where v' is the amplitude of the sound particle velocity field. 
The corresponding sound pressure field becomes. 


^ = -P A 
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Since the pressure is transmitted through the boundary- layer with- 
out change (this is one of the major benefits of boundary- layer 
theory) , Eq. (49) may be rewritten as 

The solution to Eq. (52) that satisfied the boundary condition 
given by Eq. (50) and the no-slip wall boundary conditions, v^^' = 0 
at y=0, is 
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An estimate of the fluctuating drag D' follows from Eq. (53) upon 
substitution into Eq. (47) to yield 

D'~/i ^^^<5 = + (54) 

Comparing Eq. (54) and (47), the characteristic thickness of 

the fluctuating shear stress is 
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For Eq. (55) to be valid, the characteristic length Lc}j<<6, the 
cylinder diameter. This restricts the frequency f to the follow- 
ing regime. 
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Substituting the maximum values of 6 = 10"“^ meter (Specimens #1, 5, 

6, 7) from Table IV into Eq. (56), f must exceed the minimum 
value (for the ten specimens) of f>>240 Hz. Since the frequencies 
at which specimens #1, 5, 6 and 7 were tested never exceeded 
f=1200 Hz, the effects of frequency should be small in accord with 
the measurements shown in Figure 26. Thus in light of the above 
analysis, the agreement between a.c. and d.c. resistances is to be 
expected. At very high frequencies, however, the a.c. resistance 
should become much larger than the d.c. resistance. 

5. CONCLUSIONS 

The results of this study have immediate application to the 
control of internally generated noise within jet engines. Steady- 
state and acoustic tests were performed on seven Feltmetal and 
three Rigimesh porous lining materials. The principal findings 
of this study are summarized in terms of the material steady- state 
and acoustical behavior. 
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5 . 1 Steady-State Flow Resistance 


(1) The effects of grazing flow on the resistance of thin 
porous -faced liners were measured and found to be generally 
small except for specimen #7. 

(2) The data for the seven Feltmetal specimens were corre- 

lated in terms of a semi-empirical fluid mechanical model that 
treated the porous material as a collection of interacting small 
cylinders. The correlation, shown below, is not general but 
does indicate the qualitative behavior of the resistance in 
terms of the material and medium properties. 


(3) The nonlinear behavior of the specimens was correlat- 
ed in terms of a fiber diameter based Reynolds number Re=Vp5/Pv. 
In general, the onset of nonlinearity occurred when Re ranges 
from 2 to 3. Both the data and the correlation indicate that 
small diameter fibers delay nonlinearity to higher specimen 
through flow speeds. 

5 . 2 Acoustic Impedance 

(1) The effects of grazing flow on the resistance and 
reactance of thin porous-faced liners were measured and found to 
be generally small. 

(2) A simple lumped element model was derived which 
correlated the face-sheet reactance of the ten specimens. The 
correlation, shown below, was derived by treating the porous 
material as an equivalent orifice of area PA where A is the 
cavity backing area and P is the material porosity, 

Xf o? pu) (r - 0.52 o)y/ p 


A simple prediction expression for the impedance of the seven 
Feltmetal specimens was derived based on the substitution for 
sufficiently low frequencies of the a.c. resistance by the d.c. 
resistance. An analysis presented in Section 4.3 showed that the 
a.c. and d.c. resistances are approximately equivalent for 
frequencies f<<v/2Tr6^ where v is the kinematic viscosity of the 
fluid and 6 is the mean diameter of the Feltmetal material. The 
resulting expression is 
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p 


An important finding of this study is that a great deal can 
be learned about the acoustic properties of porous lining mater- 
ials in a jet engine duct environment by the much simpler study 
of their nongrazing flow, steady-state flow resistance behavior. 
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TABLE I. SUMMARY OF JMNUFACTURERS PREDICTED ACOUSTIC PROPERTIES 


Spec. Trade Name (Nominal) NLF* Average Porosity Average 


# 

or 

Identifica- 

tion 

Predicted 
Resistance 
(Rayls -mks) 

500/20 

Thickness 

(meters) 


Pore 

Diameter 

(meters) 

Feltmetal Materials 







1 

FM 134 

320@0 . 2m/ sec 

4.7 

8.9x10''* 

0.59 

5.9x10'5 

2 

347-10-20- 

AC3A-A 

100 

2.0 

4.1 

M 

0.97 

4.8 

tt 

3 

12581-4 

580 " 

NA** 

9.1 

If 

0.80 

3.2 

tl 

4 

FM-122 

500 " 

1.8 

7.6 

M 

0.91 

3.9 

M 

5 

FM-134 

320 

4.7 

8.9 

M 

0.59 

5.9 

ft 

6 

FM-134 

320 •' 

f f 

U 

M 

ft 

It 

tl 

7 

FM-134 

320 

11 

11 

M 

II 

II 

If 

Rigimesh Materials 








8 

A0300S0204L 

2000 . 6m/sec 

NA 

3.8 

ft 

NA 

1.8 

II 

9 

A0300S0304L 

300 

NA 

3.8 

tf 

NA 

II 

If 

10 

A0300S0404L 

400 

NA 

3.8 

ft 

NA 

It 

tt 

* 

Material Nonlinear Factor 
tances at flow velocities 

defined as the ratio 
of 5.0 m/sec and 0.2 

of flow 
m/sec . 

resis 


Information not available 
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Table II. EMPIRICAL FIT OF D.C. FLOW RESISTANCE DATA 

[R = a+bVo (Rayls-mks) ] 


Specimen # 


V =0 

oo 


Blowing 


Suction 


X,= 79.6 m/s 

CO ' 


V =0 


V =79.6 m/s 


1 

398 + 240Vj3 

394+250Vp 

408 + 220VJT) 

419+210VD 

2 

117+IIVd 

120+126VD 

121+lOVp 

113+lOVj) 

3 

865+IIOVd 

860+120Vp 

855+lOOVp 

847+90VD 

4 

598+65Vd 

581+67Vd 

595+54VD 

567+61VD 

5 

493+330Vp 

489+330Vp 

505+290Vp 

563+250VD 

6 

540+320VD 

538+320Vp 

564+340Vp 

570+330Vp 

7 

637+500Vo 

644+500Vp 

717+430Vd 

733+420VD 

8 

292+36Vd 

297 + 37Vj) 

316+32VD 

298+29Vd 

9 

558+64Vp 

524+76VD 

584+61VD 

572+48Vp 

10 

298 + 34V£) 

296+35VD 

313 + 28Vj) 

282+29VD 
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TABLE III. SUMMARY OF CORRECTION TO AP_ 
GRAZING FLOW 


Specimen # 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


Corrected APg ( 
V =56.39 meters/sec 

CO 

(mks) 

.0011 m-HzO 
. 0020 
. 0020 
. 0020 

negligible 
. 0013 
. 0013 
.0018 

negligible 
. 0005 


MEASUREMENTS DUE TO 


-HzO) at Vd=0 

V =79.55 meters/sec 
(mks) 

.0023 m-HzO 
. 0038 
. 0041 
. 0043 

negligible 
. 0025 
. 0025 
.0036 

negligible 
. 0010 
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TABLE IV. CORRELATION OF D.C. FLOW RESISTANCE MEASUREMENTS 


Spec . 
# 

Rd. c 
(Rayls 
mks) 

Porosity 

- CP) 

Fiber Dia. 
(meters) 

Corr. Para- 
meter* (3) 

Obs. Vp 
at Onset 
of NL 
(meters/ 
sec) 

Values of 
Vp where 
Re=2;3 
(meters/ 
sec) 

1 

398 

. 51 

10 

1784 

.3- .4 

.15; .23 

2 

117 

.97 

8x10 

2026 

3-4 

3.63;5.44 

3 

865 

. 80 

8x10 

1538 

.7-1.1 

3 ; 4 . 5 

4 

598 

.92 

4x10 

1737 

2 . 5-3. 8 

6.92;10.4 

5 

493 

. 48 

10 

1915 

. 35- .40 

.15; .23 

6 

540 

.46 

10 

1907 

. 25- . 35 

.16; .24 

7 

637 

.45 

10 

2145 

.15-. 40 

.16; .25 

8 

292 

-.37 

4x10 

- 

.7-2 

.47; .71 

9 

558 

. 29 

4x10 

- 

1-1.5 

. 53 ; . 80 

10 

298 

.43 

4x10 

- 

.7-2 

.43; .64 

*The 

values 

of 3 tabulated above 

are determined 

from Eq. 

(41) 


with a=0 . 8 
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TABLE V. COMPARISON OF NON-GRAZING FLOW STEADY STATE, TWO MICROPHONE 
AND STANDING WAVE TUBE MEASUREMENTS 


Face Sheet Reactance 

Specific Resistance (Rayls-mks) (Xf/pwD) 


Specimen # Steady Flow Two-Mic* 


1 

400 

405 

2 

115 

141 

3 

890 

931 

4 

600 

684 

5 

500 

515 

6 

540 

553 

7 

650 

676 

8 

280 

321 

9 

560 

585 

10 

280 

235 


Impedance Tube Two-Mic* Impedance Tube 


395 

.80 

.72 

130 

.59 

.28 

906 

.69 

.40 

710 

.54 

.77 

548 

1.14 

.81 

512 

1.13 

.79 

687 

1.32 

.74 

359 

.79 

1.16 

528 

.98 

.90 

285 

.64 

.68 


*data corrected by method of Appendix A 



TABLE VI. EXPERIMENTAL DETERMINATION OF INERTIAL CONSTANT e 


Specimen # 

Thickness 
T (m) 

Porosity 

P 


'Xf 

pwD 

meas . 

e 

1 

9.6x10 

.51 


80 

.55 

2 

9.6 " 

.97 


59 

.56 

3 

9.1 " 

.80 


69 

.59 

4 

7.8 ” 

.92 


54 

.50 

5 

9.2 ” 

.48 

1. 

14 

.76 

6 

9.1 " 

.46 

1. 

13 

.73 

7 

9.1 " 

.45 

1 . 

32 

.86 

8 

1.3 " 

.37 


79 

.48 

9 

1.1 ” 

. 29 


98 

. 52 

10 

1.4 " 

.43 


64 

.43 
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STAINLESS STEEL TYPE 430 



Surface appearance at X 24 of 430 stainless steel FELTMETAL fiber 
metal. The material is 15% dense and employs Type B46 fibers. 


Photographs of Feltmetal Porous Materials 


FIGURE 1. ENLARGED PHOTOGRAPHS OF 
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Enlarged Section 
x5 Magnification 


Photograph 


of Rigimesh Porous Material 



IMEN INTERIOR AND MODEL GEOMETRY 





20 H.P. Motor 



FIGURE 2. SCHEMATIC OF STEADY-STATE FLOW RESISTANCE MEASUREMENT SYSTEM 
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FIGURE 3. TEST SECTION BOUNDARY LAYER VELOCTIY PROFILE 
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FIGURE 4. EFFECT OF GRAZING FLOW ON THE STEADY-! 
1-4 OPERATING IN THE BLOWING MODE 
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FIGURE 5. EFFECT OF GRAZING FLOW ON THE STEADY- STATE RESISTANCE OF 
SPECIMENS 1-4 OPERATING IN THE SUCTION MODE 
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FIGURE 9. EFFECT OF GRAZING FLOW ON STEADY-STATE RESISTANCE OF SPECIMENS 
8-10 OPERATING IN THE SUCTION MODE 
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10. SCHEMATIC OF VORTICAL VELOCITY INDUCED ACROSS 
SPECIMEN BY GRAZING FLOW SHEAR STRESSES 
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FIGURE 11. EFFECT OF GRAZING FLOW ON STEADY-STATE RESISTANCE OF SPECIMEN-BLOWING 





FIGURE 12. SCHEMATIC OF TWO-MICROPHONE MEASUREMENT SYSTEM 
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FIGURE 13. TWO-MICROPHONE METHOD INSTRUMENTATION 
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FIGURE 14. EFFECTS OF GRAZING FLOW AND SOUND PRESSURE LEVEL 
ON THE ACOUSTIC RESISTANCE OF SPECIMENS #1-4 
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16. EFFECTS OF GRAZING FLOW AND SOUND PRESSURE 
LEVEL ON THE ACOUSTIC RESSITANCE OF 
SPECIMENS #8-10 
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FIGURE 17. EFFECT OF GRAZING FLOW AND SOUND PRESSURE 
LEVEL ON FACE-SHEET REACTANCE OF SPECIMENS 
#1-4 
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Phase Change Across Specimen Face-Sheet 
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FIGURE 19. EFFECT OF GRAZING FLOW AND SOUND PRESSURE 

LEVEL ON PHASE ANGLE CHANGE ACROSS THE FACE- 
SHEET OF SPECIMENS #1-4 
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Non-dimensionalized Face-Sheet Reactance (X£)/pa)D 
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FIGURE 20. EFFECT OF GRAZING FLOW AND SOUND PRESSURE 
LEVEL ON FACE- SHEET REACTANCE OF SPECIMENS 
#5-7 
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Acoustic Reactance - Rayls(mks) 
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FIGURE 21. EFFECT OF GRAZING FLOW AND SOUND PRESSURE 
LEVEL ON THE TOTAL ACOUSTIC REACTANCE OF 
SPECIMENS #5-7 
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Phase Change Across Specimen Face-Sheet 
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FIGURE 22. EFFECT OF GRAZING FLOW AND SOUND PRESSURE 
LEVEL ON THE PHASE ANGLE CHANGE ACROSS THE 
FACE-SHEET OF SPECIMENS #5-7 
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Non-dimensional ized Face-Sheet Reactance - (Xfj/pmD 
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FIGURE 23. EFFECT OF GRAZING FLOW AND SOUND PRESSURE LEVEL 
ON THE FACE -SHEET REACTANCE OF SPECIMENS #8-10 
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Acoustic Reactance - Rayls(mks) 
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FIGURE 24. EFFECT OF GRAZING FLOW AND SOUND PRESSURE 
LEVEL ON THE TOTAL REACTANCE OF SPECIMENS 
# 8-10 
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Phase Change Across Specimen Face-Sheet 
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FIGURE 25. EFFECT OF GRAZING FLOW AND SOUND PRESSURE 
LEVEL ON PHASE ANGLE CHANGE ACROSS THE 
FACE-SHEET OF SPECIMENS #8-10 
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FIGURE 28 


EFFECT OF VARIATION OF FREQUENCY ON THE 
ACOUSTIC RESISTANCE OF SPECIMEN #7 
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FIGURE 29. EFFECT OF VARIATION OF FREQUENCY ON THE 
ACOUSTIC REACTANCE OF SPECIMEN #1 
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FIGURE 31. EFFECT OF VARIATION OF FREQUENCY ON THE 
REACTANCE OF SPECIMEN #7 
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FIGURE 32. EFFECT OF FREQUENCY VARIATION ON THE 
PRESSURE RATIOS AND PHASE CHANGES 
ACROSS SPECIMENS 1, 2 5 7 



APPENDIX A 


CORRECTION TO MEASURED INCIDENT SOUND PRESSURE 
LEVEL DUE TO ACOUSTIC RADIATION FROM THE TEST 
SPECIMEN 


To measure the acoustic impedance o£ a sheet o£ porous 
material backed by a cavity using the two-microphone method, it 
is necessary to accurately determine the di££erence in sound 
pressure level and phase between a microphone mounted £lush with 
the sur£ace o£ the test specimen and another mounted at the back 
o£ the resonator cavity. As has been pointed out in Section 
4.1, it is necessary to mount the "incident" microphone su££i- 
ciently £ar £rom the test specimen so that it is not in£luenced 
by the hydrodynamic near £ield. 

The typical behavior o£ highly resistive specimens was 
as £ollows. The £requency o£ the incident sound £ield was ad- 
justed to resonance which corresponds to a relative phase across 
the specimen c|) = 90 degrees. During this adjustment the incident 
sound pressure level was set at 90 dB ; was observed to de- 
crease monitonically with increasing sound pressure level. This 
was the usual response o£ highly resistive porous specimens. 
However, £or particularly low resistance specimens, c() was ob- 
served to initially increase by several degrees with increasing 
sound pressure level be£ore decreasing. Additionally, it was 
noted that in the range o£ this anomalous phase behavior, in- 
creasing the drive level to the loudspeaker by 5 dB would re- 
sult in somewhat greater than 5 dB increase in the measured inci- 
dent sound pressure level. 


This can be explained in terms o£ acoustic radiation £rom 
the test sample as £ollows: Let Q be the inward volume velocity 

in the ori£ice or porous sample, Pcav the cavity sound pres- 
sure, Pine be the incident sound pressure, (Zp) = (Rp) + (Xp) 
be the speci£ic acoustic impedance o£ the sample, normalized to 
the cavity diameter D, L be the cavity depth, and r be the 
distance £rom the acoustic center o£ the specimen to the "inci- 
dent" measuring location. 


For sinusoidal excitation Pine ~ Pinc®^^^’ where hence- 
£orth the ei^^t will be assumed. We therefore have, £or the 
volume velocity 


Q = fine 

Zd 


7TD^ 

4 


(A-1) 
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and from mass conservation and adiabatic compression. 


4pc^ 
D^L 


Qdt = -i 


4pc^Q 
7T D* LO) 


CA-2) 


Combining Eqs. (A-1) and (A-2), 


Z D L 


(^D = 


(A-3) 


The radiated pressure Pj. for a small source in a plane surface is 
given by 


Pr(^r)= -i 


pCJQ 
2 TT r 


e 


(jjr 

c 


(A-4) 


Combining Eqs. (A-1) and 
Pr(r) = -i 

Zd- 


(A-4) 


-i tJ r/c 
e ' 


Pr (r) _ pu) ^ _ i 4- 7T/ ^ 

Pine 8ZQ.r 


CA-5) 


Therefore, the total sound pressure P-t at the incident microphone 
will be given by 


Pt = P 


I nc 



1 + 


^2 

PW D 

87(r^T^ 



- I 



)1 


(A-6) 


Unfortunately, this otherwise convenient expression includes Rp 
and Xp, which we are attempting to measure, and an iterative 
calculation is required to correct for the effect. Transposing 
we have 
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-1 


^inc = 


1 + 


pu)0‘ 


8r(RD^+ Xd^)V2 



7T , Xd\ 
-;5-+arctan-;^ ) 

p*Dy 


CA-7) 


or 


Pine = A^+2 Acos^ 


i arc tan 


AsinB 
1 + AcosB 


(A-8) 


where 


A = 


ptJD^ 


1/2 


(A-9) 


and 


B = iiil- + ~ + arc tan 
c 2 



(A-10) 


Thus , 


SPLj = SPLt -10 log(j + A^ + 2 A cos B^ 


0 j = 0^ + arc tan 


( AsinB \ 
\1 + A cos B / 


(A-11) 


The iterative calculation proceeds as follows; 

1. Use SPL-fSPLi and (fif 4>i to determine on initial value of Rp 
and Xp. 

2. Use these values of Rp and Xp in Eqs. (A-8) thru (A-11) to 
determine a first iteration value of SPLj^ and cj). 
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3. Calculate new values of Rp an<J Xj). 

4. Repeat steps 2 and 3 until the answers from two successful 
iterations are equal within the desired accuracy (0.1% re- 
quires only 3-4 iterations). 

Results of this correction procedure for low sound pressure 
level data from porous samples 1 through 10 are shown in Table 
A-1. Note that for each case, ASPL is lower than measured, 

(j) is higher than measured, R/pc is slightly higher than measured 
and X/pc is more positive than measured. 
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TABLE A-1 SUMMARY OF CORRECTED IMPEDANCES TO TEN POROUS SPECIMENS 


= 0; SPL(inc) = 90 dB 


Sample 

# 

Meas 

ASPL 

CdB) 

Meas 

'f'ic 

(deg) 

Init 

R/pc 

Init 

X/pc 

Corr 

ASPL 

CdB) 

Corr 

'f'ic 

(deg) 

Corr 

R/pc 

Corr 

X/pc 

1 

.58 

90 . 1 

.9794 

.0017 

.37 

92 . 3 

1.0026 

.0397 

2 

10.07 

90.0 

. 3219 

. 0000 

9.35 

96.47 

. 3474 

.0394 

3 

-6.90 

90 . 4 

2.2791 

.0159 

-7.00 

91.37 

2.3045 

.0551 

4 

-4.01 

90.0 

1.6503 

. 0000 

-4.26 

91.53 

1.6983 

. 0455 

5 

- .82 

89.9 

1 . 2453 

- .0022 

-1.03 

91 .75 

1.2748 

.0389 

6 

-1 . 84 

89.9 

1.3400 

- . 0023 

-1.97 

91.41 

1.3595 

. 0335 

7 

-3.02 

89 . 8 

1.6500 

- . 0058 

-3.16 

91.17 

1.6770 

.0342 

8 

2 . 83 

90.1 

. 7773 

. 0014 

2.60 

92.68 

. 7971 

. 0373 

9 

7.63 

90.0 

.4725 

.0000 

7.01 

94.82 

. 5059 

. 0427 

10 

5 . 56 

89.7 

.5621 

- . 0029 

5.23 

93.30 

. 5829 

. 0336 
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